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Abstract 

Two basic classes of electromagnetic media, recently defined and 
labeled as those of P media and Q media, are generalized to define 
the class of PQ media. Plane wave propagation in the general PQ 
medium is studied and the quartic dispersion equation is derived in 
analytic form applying four-dimensional dyadic formalism. The result 
is verified by considering various special cases of PQ media for which 
the dispersion equation is known to decompose to two quadratic equa¬ 
tions or be identically satisfied (media with no dispersion equation). 

As a numerical example, the dispersion surface of a PQ medium with 
non-decomposable dispersion equation is considered. 

1 Introduction 

The most general linear (bi-anisotropic) electromagnetic medium can be rep¬ 
resented in terms of 36 scalar medium parameters in different representations. 
In engineering form applying Gibbsian vector fields and 3D dyadics it is typ¬ 
ical to apply the form [1, 2] 
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Another form favored by the physicists is [3] 



A more compact form is obtained by applying the 4D formalism in terms of 
differential forms [4, 5] where the electromagnetic helds are characterized by 
two six-dimensional quantities, the two-forms and 

= D - H A £ 4 , # = B -f E A £ 4 . (3) 

In a linear medium they are related by the medium bidyadic M as 

= (4) 

The medium bidyadic can be expanded in terms of 3D medium dyadics as 

M = tt -|- 6 A 64 -|- £4 A yU -|- £4 A /3 A © 4 . (5) 

For details in the present notation [5] or [ 6 ] should be consulted. 

It is not easy to get a grasp of the most general medium dehned by the 
bidyadic M. This is why many classes of media with M restricted by special 
forms involving less than 36 parameters, and bearing strange names, have 
been defined and studied. Because M corresponds to a 6 x 6 matrix, some 
classes have been based on expressing the bidyadic M in terms of a dyadic 
corresponding to a 4 x 4 matrix involving only 16 parameters. As the most 
obvious one, the class of P media [7] has been dehned in terms of a dyadic 
P G EiFi as the double-wedge square 

M = (6) 

As another obvious example, the class of Q media [ 8 ] has been dehned 
through the modihed medium bidyadic 

Mm = eAr[M = eAr[Q*'^^ G E 2 E 2 , (7) 

in terms of a (quasi-) metric dyadic Q G EiEi. Properties of both medium 
classes have been recently studied. For a plane wave, a held of the depen¬ 
dence exp(n'|x) of the spacetime variable x, the wave one-form jz is normally 
restricted by a dispersion equation 

D{jy) = 0, 


( 8 ) 
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which is a quartic equation in general [9] and its form depends on the medium 

bidyadic M. For example, for the Q medium the dispersion equation is of the 
form 

D{v) = ^Q{Q\\uuf = (9) 

which means that, for a dyadic of full rank satisfying Ag = 7 ^ 0 , 

the quartic equation reduces to a quadratic equation. On the other hand, 
one can show that, for a P medium, the dispersion equation ( 8 ) is actually an 
identity which is satisfied by any one-form v [7]. Also, if Q is of rank lower 
than 4, we have Ag = 0 and (9) is satisfied identically. Media of this kind 
have previously been called NDE media (media with no dispersion equation) 
| 10 |. 

As simple generalizations of the classes of P and Q media we can define 

M = (10) 

and 

M^ = W^^ + AB, ( 11 ) 

which have respectively been called extended (or generalized) P [7] and Q 
[11] media. Here, A, B, C and D are any bivectors. The nnmber of medium 
parameters in both cases has been increased from 16 to 23. 

Plane waves propagating in extended P and Q media have been previously 
studied and the dispersion eqnations have been shown to take the respective 
form [7, 11] 


D{i.) = Ap{(p[P^)\^u){{C[P-^^)\\uu) =0, ( 12 ) 

D{iy) = Ag(Q||jzjy)(QAB[[Q~^'^)||jzjy = 0. (13) 

Here we denote Ap = trP’^'^b In both cases, for full-rank dyadics P, Q, the 
dispersion equations can be decomposed in two quadratic equations. Such 
media make two examples of what have been called decomposable media [ 12 ]. 
For dyadics P and Q of rank less than 3, the two media again fall in the class 
of NDE media with dispersion equations satisfied identically for any jz [10]. 
For ranks equaling 3, some qnadratic fnnctions of D{iy) may decompose to 
prodncts of linear functions of u. 
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2 Class of PQ Media 

In the present stndy we make a further generalization of both P and Q media 
by considering medium bidyadics of the composite form 

W = + Sn[W^^. (14) 

Media dehned by bidyadics of this form will be called PQ media. It is easy to 
see that medium dyadics of this form do not cover all possible linear media. 
In fact, the number of parameters must certainly be less than 2 x 16 = 32, 
which falls short of the 36 parameters of the most general medium. 

To hnd a 3D representation of the PQ medium bidyadic (14) the dyadics 
P and Q can be expanded in terms of spatial vectors as,b 5 ,p^, a spatial 
one-form TZg and spatial dyadics Q^, Ps as 


Q — Qs -l- G4as -|- bse 4 -|- 06464 , (15) 

P = P 5 -h 6477 ^ -h p ^£4 + ^ 6464 . (16) 

The spatial medium dyadics of the Q medium can be expressed as [5] 

^ = Si23LQsAa3, (17) 

^ = Si23L(cQs - b^aQ, (18) 

= -^123 (19) 

/3 = ^123L(bs A Qs), (20) 

while those of the P medium have the form [7] 

S _ (21) 

r = -^.APf, (22) 

= -Pf A Ps, (23) 

^ = TTsPs-pPr. (24) 


The medium dyadics of the PQ medium are obtained by summing up the 
corresponding expressions. To write the medium equations in the Gibbsian 
form (1) requires that the dyadic J1 exist. Denoting the Gibbsian dyadics by 
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the subscript ()g, they are obtained by inserting the previous expansions in 
the expressions [5] 



= ei23L(e'- 

(25) 


= 0123 [a /i) 

(26) 

Cg 

= —ei23L/^F) 

(27) 

f^g 

= ^123 [fj,. 

(28) 


The dyadic denoted here by Jl stands for the inverse of the sum of the two 
dyadics of (19) and (23). The full analytic form of the Gibbsian medium 
dyadics of the PQ medium would have quite an extensive form. 

3 Plane Wave in PQ medium 

The main task of this study is to hnd properties of a plane wave propagating 
in the general PQ medium. Expressing the held two-form of a plane wave in 
terms of a potential one-form as ^ = u A (p, the potential satishes 

u A ^ = u A MKu A <p) = 0 

= uA{W^^^ + eN[W^'>)\{iyAcl>) 

= A (P^|j^) A (P^|0) -h SAr[(Q^^^[[j/J/)|0. (29) 

Operating this by eAr[ yields the equation 

e7vL(i^A^') = A{iy\P) AP'^)\<p+ 

Here we have introduced the bivector 

F = F(z.) = e^L(^A(z.|P)), (31) 

which is simple since it satishes [6] 

FAF = 0, (F[F)(2) = FF, (F[F)(3)=0, (32) 

and it can be expressed in the form F = a A b in terms of two vectors. The 
equation (30) for the potential, 


D{v)\4i = 0 , 


( 33 ) 
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is defined by the dispersion dyadic 

Because of (33) and D(i/)|t/ = 0, the rank of D(i/) must be less than three, 
provided (p and i' are linearly independent, i.e., when # = i/ A 0 7 ^ 0, which 
is assumed here. As a consequence, u is restricted by the dyadic dispersion 
equation [ 6 ] 

5(3) (jy) = + 

= (F[F)(3)|P(3)^ + ((F[F)(2)|p(2)r)0(5(2)LLiyiy) 

+ (FLF)0(5(2)^^^^)(2) ^ (5(2)lL^^)(3) 

= 0. (35) 

Applying the expansion rules [ 6 ] 

= (Q||j^z^)Q - (Q|z^)(z/|Q) 

= (36) 

= (Q||j/z/)^(Q(^([[j/j/) = AQ(Q||jyjy)^(ejve 7 v[[z/z/), (37) 
with Aq = £4641 |Q("^), the dispersion equation (35) can be written as 


D(3)(iy) = Ci(z/) + € 2 ( 1 ^) + Csi^jy) = 0, (38) 

where we denote 

Fi(F = (FF|F2)'^);((5^2)^^^^) (39) 

F 2 (F = (^||*^F(fLF^)0(^®LL^F (40) 

€ 3 ( 1 ^) = AQ{Q\\uuf{eNeN[[iyiy)- (41) 


Now one can show that (38) is equivalent to a scalar dispersion equation ( 8 ). 
For that we expand the three dyadics Cpiy) as follows. 
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• For the dyadic Ci(i/) we apply the identity 

AiA(Bi[Q:) = (Aj A Bi) [cK - Bi A (A* [ck) 

= £N\iAiABi){eN[a)-BiA{Ai[cx), (42) 

valid for any bivectors Aj,Bj and one-form a. Assuming that Aj and 
(X satisfy Aj[Q: = 0 , we can construct the dyadic rule 

AiA 2 a(BiB 2 [[cKCk) = £Ar£Ar||(AiA2ABiB2)(ejve7v[[Q:Q:)- (43) 

Because of F[i/ = 0 and (F|P(^)'^) [i/ = (F[(i/|P)) = 0, we can set 

Ai = F, A 2 = F|P*^^^^ and cx = u. Since the rule is linear in the dyadic 
B 1 B 2 , we can set B 1 B 2 = and apply the rule as 

Zi{u) = (FF|F(2)^)0(Q(2)[[j/j/) = Di{u){e]^eN[[jxJx), (44) 

with 

= AQ^(-2)^||(FFp)^) = AqFF||(F^|^-1)(2). (45) 

In the last expression we have assumed that Q is of full rank, Aq 7 ^ 0. 

• To expand the dyadic € 2 ( 1 ^) we proceed as 

= AQ(Q||j/j/)((eAreAr[[(Q"^^0j/jy))0(F[P^) 

= Aq(Q| |z/z/)eAreAr [[((Q"^^0z/j/) [[(F [P^)) 

= AQ(Q||j/j/)(e7veAr[[jyjy)(Q"^'^||(F[F^)) 

= AQ(Q||z/z/)(e7veAr[[jyjy)(Q"^'^||(F[F^)) 

= D2{u){eNeN[[i^i^), (46) 

with 

D2iu) = Ag(^||zyiy)tr(FLF|^-i). (47) 



Finally, we have 


(48) 


Ci(v) = D3(i/)(e^,eN[[i/i/), 

with _ 

Dsiiy) = AgiQUuiyf. (49) 

Because each of the dyadics Q(i/) is a scalar multiple of otvStv the 
dyadic dispersion equation (38) equals the scalar dispersion equation (8) as 

D{u) = Di{u) + D2 {u) + D3{u) 

= AqFF||(?^|^- 1)(2) + AQ(^||iyi/)tr(FLP^|^-') 

+AQ(^||iyiy)2 = 0. (50) 

Substituting F = eAr[(i/A (r/|P)), the quartic form of the dispersion equation 
can be . This is the main result of the present paper. 

For Aq —)■ 0 we must replace 

AqQ*- AqQ ^ > £Ar£Ar (51) 

in the expression (50). 

4 Special Cases 

Let us consider the expression (50) for a few special cases of the PQ medium 
for which we know the dispersion equation. 

1. For the pure P medium case Q = 0, after inserting (51) we obtain the 
identity D{iy) = 0 for all u. This proves that a pure P medium belongs 
to the class of NDE media [10]. 

2. For the pure Q medium with P = 0 (50) is reduced to the known 
quadratic dispersion equation (9). 

3. The case P = al corresponds to a Q-axion medium. Since we now have 

F = ae]\f[{h' A u) = 0, only the last term of the expression in (50) 
survives. This, again, yields the dispersion equation of the Q medium. 
In fact, it is well known that adding an axion term to the medium 

bidyadic does not change the dispersion equation. Here one should note 
that the P-axion medium M = + with a ^ 0 is not a special 

case of the PQ medium. 
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4. Choosing Q = aibi + a 2 b 2 we obtain 

M = p(2)'r g.^|^(ai A a2)(bi A b2), (52) 

which yields a special case of the extended P medium (10) with DC = 
(ai Aa 2 )(bi Ab 2 ). Applying (51) for Aq —)■ 0 in (50) yields the disper¬ 
sion equation whose form can be expanded as 

D{u) = 

= F|F(2)^|(s^L(bi A b2))((ai A a2)JsA,)|F 
= {i/ A (j/|P))Je7v|P^^^'^|£iv[(bi A b2))((ai A a2)|(i^ A (i^|P)) 
= Ap{u A (j/|P)|p("^^|(bi A b2))((ai A a2)|(j^ A (j^|P)) 

= Ap((jy|?“^) A J^)|(bi A b2))((ai A a2)|(j^ A (j^|P)) = 0.(53) 

Assuming Ap = trP^^^ ^ 0 the equation is split in two quadratic equa¬ 
tions as 

j^|(P A j^)|(ai A as) = j/|(PJ (ai A a2))|j/ = 0, (54) 

A J^)|(bi A bs) = J/|(P'^J(bi A bs))]!^ = 0, (55) 

which coincide with (12) for DC = (ai A a 2 )(bi A bs). In the case 
Ap = 0 we must replace P~^ by eArSAr in (55). 

5. Choosing P = aicci -|- asccs we have P*-^^ = 0 and 

M = SatLQ^^^^ + («! A Q:2)(ai A as), (56) 

which corresponds to a special case of the extended Q medium (11) 
with AB = eAr[(Q:i A Q: 2 )(ai A as). Expanding 

F|P(2)T = A = ej^liu A = 0, (57) 

the hrst term of (50) obviously vanishes. Expanding further 

F[W = e^L(^A(F|i/)AF) = e^L(^A(p2)^L^)) 

= -u\ {bn [P^^^^) [jy = {bn [P^^^^) [ [I'l', 


( 58 ) 
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we have 




whence the dispersion equation (50) can be written as 


D{u) = /\q{Q\\vv){Q + (e^[p(2)^)[[Q- 
which coincides with (13). 


\vi' = 0, (60) 


6 . Finally, let us assume that Q is an antisymmetric dyadic, which can be 
expressed in terms of some bivector A in the form 

^ = A[F. (61) 

This implies QWw = 0, whence the dispersion equation (50) reduces 
to _ _ 

D{v) = AqF|F2)^|Q(-2)|F = 0. (62) 

Applying the expansion [ 6 ] 

(ALF)(2) = AA + aejvJ^')^, a =-^Sjv|(A A A), (63) 

the PQ medium bidyadic has the form of an extended P-axion medium 
bidyadic [ 12 ] 

M = [AA + (64) 

Applying the rule [ 6 ] 

= ^(£Ar£Ar[[AA + a£Ar[F)), (65) 

one can show that its last term has no effect on the dispersion equation. 
In fact, inserting 

F|F2)^ = e7v|(i^A(F|r^) aF^)'^) =zy|(e^[F3)^)[zy 

= Apu\{P~'^\eN)[u, ( 66 ) 

in (62) we have 

F|P(^)'^|(£Ar[F^)|F = Ap(z^|(P”^JeAr)[j^)|(£Ar[(ejv[(j/A (P^F))) 

= Ap(r/|(P“^JeAr)|(i/A r/A (P'^F)) = 0. (67) 



11 


Thus, (62) is reduced to 

D{u) = 

= ApJ/|(P^^JeAr)[j^|(£Ar£Ar[[AA)|(e7v[(j^ A (P^|j^)) 

= ApJ/|(P^^JeAr)[j^|(£Ar[A)(A|(j/ A (P^|j^))) 

= Ap((zy|p-1 Ar^)|A)(A|(jyA(F^|jy))) 

= Ap((p-iJA)||r.r.)((FjA)||^r.) = 0. (68) 

Since the axion component does not affect the dispersion equation, the 
result coincides with that of the extended P medium (12) for the special 
case C = D = A. 

5 Example 

As a numerical example of a PQ medium let us consider a special one by 
restricting the dyadic P as 

P = P(l+Po), Po = acK, a|Q: = 0, (69) 

and assuming that Q be a symmetric dyadic, 

^ = 1, = f. (70) 

In the corresponding medium bidyadic 

M = p2|(2)T ^ p2^^Ayr ^ y|(2) 

the three terms can actually be recognized as components of the Hehl- 
Obukhov decomposition, respectively called as the axion, skewon and prin¬ 
cipal components [3, 6]. In fact, the axion part is a multiple of the unit 
bidyadic while the skewon and principal parts are trace free. Also, any 
skewon bidyadic is known to be of the form (Bq))!)^ with a skewon-free dyadic 
Bo while any bidyadic C G F2E2 satisfying C[[l = 0 is a principal bidyadic. 
Actually, one can show that (^tv LU vanishes for any symmetric dyadic 
§, [ 6 ]. 
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To find the dispersion equation for the present PQ medium from (50), let 
us hrst expand 

F|F(2)t _ e7v|(jy A 

= P^ejvK*^ A (F^ + aa;(p)'^))[zy 
= P^e^Ki/ Au A + P^e^vKi^ A (Q:(a|^') A 
+i/ A (cKa) A F) 

= P^{a\u)eN[{u Acx), (72) 

FLF = eArL(i^A(F2)^))Li^) 

= P^bn [{i' A + CKaOP) [u) 

= P'^e^Ku A A F + r/ A (Q:(a|r/) A F + r/ A aa)) 

= P^(a|z/)(e7v[(i^ A q:))[F. (73) 

Since the last expression is an antisymmetric dyadic and A 5 S~^ = eArSAr 
is a symmetric dyadic or zero, we have 

AQtr(FLF^I^-i) = A5(FLF)||!-i = 0. (74) 

Thus, (50) becomes 

D{iy) = AsF\{P^\^-Y^\F + As(E\\uuf. (75) 

Finally, applying (72), we can expand 

AsF|(P^|l~^)(^)|F = P^(a|zz)(eArL(i^ A Q:))|l^^^|(e 7 v[(jz A (F + Q:a)|zz))) 
= P^(a|r/)^(r/A q:)|S‘'^^|(iz A ck)) 

= P^(a|i/)^(Q:Q:JJS‘'^^)||jzjz, 

whence the dispersion equation for the special PQ medium (71) has the form 

D{u) = P‘^(a|r/)^(Q:Q:JJS‘'^^)||iziz + A5 (S||iziz)^ 

= P‘^(a|t/)^(Q:Q:||S)(S||i/i/) + P"‘(a|i/)^(Q:|S|iz)^ 

+As(f||^^)^ = 0. 


(76) 
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Unlike for all the special cases considered above, the quartic dispersion equa¬ 
tion corresponding to the medium dehned by (71) does not necessarily de¬ 
compose in two quadratic equations. In the special case when the dyadic S 
satishes S|q: = Aa, whence we have S||q:q: = 0, the dispersion equation is 
decomposable. 

Setting S = 0 in (76) we obtain D{i') = 0 for any v, a property shared 
by all skewon-axion media. For aa = 0 we are left with As'(S||n'n')^ = 0, 
valid for the Q-axion medium. 

To be able to work on this example numerically, let us choose 



S - S(3g — 5 G 4 G 4 , 

(77) 

where denotes the metric dyadic 



Gs — GiCi -|- 6262 -|- 6363. 

(78) 

Thus, we have 

A 5 = 

(79) 

Let us further choose 

OL = 62 , a = 63. 

(80) 


Inserting these in (76), the dispersion equation becomes 


D{i') = P^vlS{S{vl + + vl) — sk"^) + 

—S^s{S{vl + 1^1 + i^l) — sk'^Y = 0 - ( 81 ) 

For a given value of k, this corresponds to a quartic dispersion surface in the 
three-dimensional space spanned by z/j. 

Let us assume that S, s and P are real and positive quantities and simplify 
the expressions by denoting 

p4 

Xi = K = ky/s, T = (82) 

with 0 < r < 1. The dispersion equation (81) then takes the form 

(xj -f a;2 + (1 — t)xI — K^){xl + xl + xl — K"^) — tx^x^ = 0. (83) 

An idea of the dispersion surface can be obtained by considering the three 
main sections. 
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• Assuming x^, = 0, (83) leads to 

x\ + X 2 — K'^ = 0, (84) 

which corresponds to a circle of unit radius. 

• Assuming X 2 = 0, (83) yields 

{xl + (1 — t)xI — K'^){xl + xl — K‘^) = 0. (85) 

This splits in two separate curves, one of which is a circle of unit radius 

xl + xl — K‘^ = 0 , ( 86 ) 

and, the other one, a quadratic curve 

xl + {l- t)xI -K^ = 0. (87) 

For r < 1 the latter dehnes an ellipse with axial ratio v^l — r. For 
r > 1 the curve is a hyperbola. 

• Finally, assuming xi = 0, (83) yields 

_ ^ 2 ) ^^2 ^ ^2 _ ^ 2 ) _ ^^ 2^2 ^ ( gg ) 

This corresponds to a curve of the fourth order, the form of which 
depends on the parameter r. For r < 1 the curve is closed and for 
r > 1 it is open. For r —)■ 0 the curve approaches a circle of unit 
radius, in which case the PQ medium approaches a Q medium. 

The cross sections Xi = 0 and X 2 = 0 are depicted for the parameter value 
r = 0.7 in Fig. 1 in terms of normalized axis parameters nui = Xi/K. It is 
seen that, for this particular PQ medium, there is no birefringence for waves 
whose wave one-form satishes e^\i' = 0. 


6 Conclusion 

A novel class of electromagnetic media, called that of PQ media, was intro¬ 
duced as a generalization of the previously studied classes of P media and Q 
media. Plane-wave propagating in the general PQ medium was studied and 
the quartic dispersion equation was derived in analytic form. The equation 
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Figure 1: Two cross sections of the quartic dispersion surface corresponding 
to 1^1 = ei|i/ = 0 and z /2 = 62 ]^' = 0. For 1^2 = 0 the cross-section reduces to 
a circle and an ellipse. For 1^3 = 0 the dispersion surface reduces to a single 
circle of unit radius. Here the parameter value r = P'^/S^s = 0.7 has been 
assumed. 
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was verified for six special cases of PQ media for which the analytic form has 
been found from previous studies. In all of these special cases the quartic 
equation either reduces to two quadratic equations or becomes an identity. 
As an example of a medium yielding a more general quartic dispersion equa¬ 
tion, another special case of the PQ medium was considered. 
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